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Abstract 

Bilayer graphene nanoribbon with zigzag edge is investigated with the tight binding model. Two stacking structures, 
a and ft, are considered. The band splitting is seen in the a structure, while the splitting in the wave number direction 
is found in the ft structure. The local density of states in the ft structure tend to avoid sites where interlayer hopping 
interactions are present. The calculation is extended to the boron-carbon-nitride systems. The qualitative properties 
persist when zigzag edge atoms are replaced with borons and nitrogens. 
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1. Introduction 

The graphite, multi-layer, and single-layer graphene 
materials have been studied intensively, since the elec- 
tric field effect has been found in atomically thin 
graphene films [1]. These materials can be regarded as 
bulk systems. On the other hand, nanographenes with 
controlled edge structures have been predicted to have 
localized states along the zigzag edges [2]. The pres- 
ence of the edge states has been observed by experi- 
ments of scanning tunneling spectroscopy [3,4]. Thus, 
the studies of the edge states are one of the interesting 
topic of the field. The recent atomic bottom-up fabri- 
cation of nanoribbons really promotes experimental and 
theoretical investigations [5]. 

In the previous paper [6], the tight binding model has 
been solved numerically, and effects of interlayer inter- 
actions have been considered. In the ft structure [Fig. 1 
(b)], the split of the energy bands is not seen compared 
with that of the a structure [Fig. 1 (a)]. The difference 
of the band split in the a and ft structures appears. Num- 
ber of states in the energy window for the positive wave 
number has been calculated for the single layer, a, and 
ft structures. The dependence of the energy reflects the 
band structures, and this will appear in quantization of 
conductance experimentally. 

In this paper, we will report the electron density dis- 
tribution in nanographene materials with zigzag edges 
including the inter-layer interactions. The calculation is 
extended to the boron-carbon-nitride systems, too. We 
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report the band structures of the single and bilayer sys- 
tems. 



2. Model 



The following tight binding model is used in the cal- 
culation: 

H = ~ t Z (c *> c ^ +h - a) 

<i,j),cr 

-h ^(cljj-r+h.c), (1) 

<U>,<r 

where c,> and di iCr are the annihilation operators of elec- 
trons at the lattice site i of the spin <x on the upper and 
lower layers, respectively. The quantity t is the hopping 
integral of n electrons between neighboring lattice sites. 
Two stacking patterns, shown in Figs. 1 (a) and (b), 
are considered. They are named as a and ft structures, 
respectively. This convention has been used in the liter- 
ature [7]. There are N z — 10 zigzag lines in upper and 
lower layers. Inclusion of the bond alternations [8,9] 
would be interesting, but it is not treated in this paper. 
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Fig. 1 . A-B stacked bilayer graphene nanoribbon with 



zigzag edges. The upper layer is shown by the solid 
lines, and the lower layer by dotted lines. In the a 
structure (a), the upper layer is shift by the bond length 
downward to the position of the lower layer. The region 
surrounded by the dashed line is the unit cell in the di- 
rection of the one dimensional direction. At the circles, 
two carbon atoms of the upper and lower layers overlap 
completely, and there is the weak hopping interaction t \ 
here. In (b), the /3 structure is shown, where the lower 
layer is shift to the right-down direction so the stacking 
pattern is different from that of the a structure. The con- 
vention of the a and /3 structures is the same used in the 
paper [7]. In the area surrounded by the dashed lines, 
the unit cell is displayed. 

3. Electron density in the real space 

In order to understand effects of the band split clearly, 
it is useful to watch the distribution of electron density 
in the real space. Here, we calculate electron density 
by changing the energy window. Figure 2 shows local 
density of states of electrons in the unit cell of the up- 
per layer. The energy is shown below each figure, with 
the number of states N s crossing at the energy for the 
region of positive wave number. The diameter of the 
circle is proportional to the electron density. The total 
electron number in the unit cell is normalized to be con- 
stant divided by N s , in order to eliminate the effects of 
enlargement of circles due to large degeneracy. In Fig. 
2, the electron density is displayed for N z = 10, so there 
are 20 atoms in the unit cell. In the single layer (a), the 
electron density is symmetric with respect to the center 
of the unit cell. At the energy E = 0.01 f, the electron 
density is large at A type sites in the upper part, while 
it is large at B type sites in the lower part. There will 
be an edge state at E = 0, but the electron density has 
already intrude into the inner part at this energy. As 
the energy increases from E = 0.1 1 to 0.7?, the electron 
density moves gradually from A type to B type sites (or, 
from B type to A type sites). After interlayer interaction 
t\ — O.lf is taken into account, the spatial symmetry is 
lost as shown for the a (b) and f3 (c) structures. The 
energy is changed from E = O.Olf to 0.75f in (b), and 
from E = O.Olf to 0.7f in (c). N s is increased by steps of 
jumps. Both in (b) and (c), the electron density is rela- 
tively larger at odd number sites from the top of the unit 
cell. It is smaller at even number sites in the inner part. 
The result looks similar at a glance. However, there is a 
difference. The interlayer interaction t\ is present at odd 
number sites of the unit cell in the upper layer in Fig 1 . 
(a). So, this interaction gives rise bonding anti-bonding 
split of energy band. The electrons favor to have large 



amplitudes at odd number sites in the inner part. On the 
other hand, the interaction t\ is located at even number 
sites in Fig. 1 (b). The electrons tend to avoid the inter- 
action 1 1 , and the removal of the degeneracy occurs in 
the direction of the wave number. In this way, the elec- 
tron density in the real space is related mutually with the 
difference in changes, seen in the band structures of the 
bilayer nanoribbons. These findings will be observed 
possibly by conduction and STM experiments. 



(a) single layer 




E = 0.01t E = 0.1t E = 0.45t E = 0.7t 
(N s = 1) (N s =1) (N s = 3) (N s = 5) 



(b) a structure 




E = 0.01t E = 0.1t E = 0.4t E = 0.5t E = 0.65t E = 0.75t 
(N s = 2) (N s = 2) (N s = 4) (N s = 6) (N s = 8) (N s = 10) 



(c) (3 structure 




E = 0.01t E = 0.1t E = 0.45t E = 0.7t 
(N s = 2) (N s = 2) (N s = 6) (N s =10) 

Fig. 2. Local density of states in the unit cell of the 
upper layer for the single layer (a), a structure (b), and /3 
structure (c). The density is proportional to the diameter 
of the circle. The energy is shown below each figure, 
with the number of states N s crossing at the energy for 
the region of positive wave number. 

4. Band structures of the BCN systems 

The calculation is extended to the boron-carbon- 
nitride (BCN) systems [10]. The stacking patterns are 
considered with the spatial inversion symmetry. The 
presence of B and N atoms is taken into account by the 
site energies Eb = +t and E^ = -f, as has been used in 
the paper [11]. 

Energy band structures of the N z = 10 systems are 
displayed in Fig. 3, for the single layer (a), a struc- 
ture (b), and [3 structure (c). The interlayer interaction 
strength is t\ = O.lf. In (b), the split of the energy bands 
is seen compared with Fig. (a). In contrast, split in the 
perpendicular direction is small in Fig. (c), except for 
the lowest unoccupied state and highest occupied state, 
which have large amplitude near edge atoms. In order 
to see split structures clearly, details near the Brillouin 
zone edge are magnified for the single layer (d), a struc- 
ture (e), and /3 structure (f). In (e) and (f), the energy 
bands of the single layer are shown by the red lines for 



comparison. We find that the qualitative properties be- 
tween the difference of the energy band splits of the a 
and j6 structures persist when zigzag edge atoms are re- 
placed with borons and nitrogens. This property is con- 
firmed by looking at the numerical data, also. 




Fig. 3. Energy band structures of the boron-carbon- 
nitride nanoribbons for the single layer (a), a structure 
(b), and /3 structure (c). Details near the Brillouin zone 
edge are magnified for the single layer (d), a structure 
(e), and /3 structure (f). In (e) and (f), the energy bands 
of the single layer are shown by the red lines for com- 
parison. 

5. Summary 

In summary, weak interlayer interactions have been 
considered for the bilayer graphene nanoribbon with 
zigzag edge. The a and f3 stacking structures have been 
considered. The band splitting is seen in the a structure, 
while the splitting in the wave number direction is found 
in the /3 structure [6]. The local density of states in the 
/3 structure tend to avoid sites where interlayer hopping 
interactions exist. The calculation has been extended to 
the boron-carbon-nitride systems. The qualitative prop- 
erties persist when zigzag edge atoms are replaced with 
borons and nitrogens. 
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